Nucleation of Ge clusters at high temperatures on Ge/Si(001) wetting layer 
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Difference in nucleation of Ge quantum dots during Ge deposition at low (<600°C) 
and high (>600°C) temperatures on the Si(OOl) surface is studied by high resolution 
scanning tunneling microscopy and in situ reflected high-energy electron diflFraction. 
Two process resulting in appearance of {105}-faceted clusters on the Ge wetting layer 
have been observed at high temperatures: Pyramids have been observed to nucleate 
via the previously described formation of strictly determined structures, resembling 
blossoms, composed by 16 dimers grouped in pairs and chains of 4 dimers on tops of 
the wetting layer M x N patches, each on top of a separate single patch, just like it 
goes on at low temperatures; an alternative process consists in faceting of shapeless 
heaps of excess Ge atoms which arise in the vicinity of strong sinks of adatoms, such 
as pits or steps. The latter process is not observed at low temperatures; it is typical 
only for the high-temperature deposition mode. 
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I. INTRODUCTION 



Issues of nucleation of {105}-faceted Ge clusters under different growth conditions on 
Si(001_) have attracted an attention of researches since the discovery of huts^ (see, e.g., 
Refs.|2|-l6| and a brief review presented in the introduction to Ref. 7). In recent articles,-"— 
we described a process of hut nucleation at low (<600°C) temperatures of Ge deposition 
consisted in formation of strictly determined structures composed by 16 dimers grouped in 
pairs and chains of 4 dimes on tops of the wetting layer M x N patches, each on top of a 
separate single patch (Fig.[l]). However, it has been unclear thus far if the same nucleation 
process takes place at high temperatures (> 600°C) or there are some different cluster for- 
mation mechanisms peculiar for this growth mode. Previous articles by different authors^*^, 
who had explored hut formation in the process of the gas-source molecular beam epitaxy 
(MBE) based on decomposition of germane (GeH4) at high temperatures on a Si surface, 
indicated that huts arose near wetting layer (WL) irregularities, such as {105}-faceted pits, 
which appeared in thick hydrogen-rich WL prior to cluster nucleation, or [100] steps. In fact, 
it has been almost commonly adopted that huts always nucleate as a result of the latter 
process regardless the epitaxy technique applied for Ge film deposition; [100] steps of much 
thinner WL, which forms in the ultra-high vacuum (UHV), have been assumed to be the 
required irregularities in the case of UHV MBE. 

Being aware of a different hut nucleation mechanism at low temperatures we have decided 
to re-examine the outset of hut array formation at high temperatures to verify or refute the 
commonly adopted opinion. The results of our investigations are presented in the current 
article. 



II. EQUIPMENT AND EXPERIMENTAL TECHNIQUES 

The experiments were carried out using an integrated UHV instrument^^^iiii^ built on 
the basis of the Riber SSC2 surface science center with the EVA 32 MBE chamber con- 
nected through a gate to the STM GPI-300 UHV scanning tunneling microscope^^. A 
preliminary annealing and outgassing chamber, as well as facility for reflected high energy 
electron diffraction (RHEED) analysis and a number of complimentary techniques of surface 
exploration, are also available in the instrument.— 
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The pressure of about 5 x 10~^ Torr was kept in the prehminary anneahng chamber. The 
MBE chamber was evacuated down to about 10~^^ Torr before processes; the pressure grew 
to nearly 2 x 10~^ Torr at most during the sample surface deoxidization process^^^^^ and 
10~^ Torr during Ge or Si deposition. The residual gas pressure did not exceed 10~^^ Torr 
in the STM chamber. 

Sources with the electron beam evaporation were used for Ge deposition. The deposition 
rate and coverage were measured using the Inficon Leybold-Heraeus XTC 751-001-Gl film 
thickness monitor equipped with the graduated in-advance quartz sensors installed in the 
MBE chamber. Tantalum radiators were used for sample heating from the rear side in both 
preliminary annealing and MBE chambers. The temperature was monitored with chromel- 
alumel and tungsten-rhenium thermocouples of the heaters in the preliminary annealing and 
MBE chambers, respectively. The thermocouples were mounted in vacuum near the rear 
side of the samples and in situ graduated beforehand against the IMPAC IS 12-Si pyrometer 
which measured the sample temperature through chamber windows. The composition of 
residual atmosphere in the MBE camber was monitored using the SRS RGA-200 residual 
gas analyzer before and during the process. 

The STM tip was ex situ made of the tungsten wire and cleaned by ion bombardment— 
in a special UHV chamber connected to the STM one. 

In this work, the images were obtained in the constant tunneling current (/t) mode at 
the room temperature. The STM tip was zero-biased while the sample was positively or 
negatively biased (Us) when scanned in empty-states or fiUed-states imaging mode. Original 
firmwareiSii^ was used for data acquisition; the STM images were processed afterward using 
the WSxM software^'^. 

RHEED analysis was performed in situ during sample annealing or Ge deposition by the 
Staib Instruments RH20 diffractometer. 

Initial samples for STM were 8x8 mm^ squares cut from the specially treated commercial 
boron-doped Czochralski-grown (CZ) Si(lOO) wafers (p-type, p = 12 J~icm). After wash- 
ing and chemical treatment following the standard procedure described elsewhere^^ (which 
included washing in ethanol, etching in the mixture of HNO3 and HF and rinsing in the 
deionized water— i^^), the silicon substrates were mounted on the molybdenum STM holders 
and infiexibly clamped with the tantalum fasteners; the STM holders were placed in the 
holders for MBE made of molybdenum with tantalum inserts. Then, the substrates were 
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outgassed at the temperature of around 565°C for more than 6 h in the prehminary anneahng 
chamber. After that, the substrates were moved into the MBE chamber where they were 
subjected to final anneahng with the maximum temperature of about 925°G^^^^. Then, 
the temperature was rapidly lowered to about 750°C. The rate of the further cooling was 
around 0.4°C/s that corresponded to the quenching mode applied in Ref.lll 



The resultant 

Si (001) surfaces were completely purified of the oxide film and atomically clesiio^^^^] their 



structure was investigated in detail in Refs. 
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14. 



Ge was deposited directly on the deoxidized Si(OOl) surface; the deposition rate {dhce/dt) 
was varied from about 0.1to0.15A/s. The substrate temperature during Ge deposition (Tgr) 
was 360°C for the low-temperature mode and 600 or 650°C for the high-temperature mode. 
The rate of the sample cooling down to the room temperature was approximately 0.4°C/s 
after the deposition. After cooling, the prepared samples were placed for analysis into the 
STM chamber. 



III. NUCLEATION OF GERMANIUM CLUSTERS AT LOW AND HIGH 
TEMPERATURES 

A. Scanning Tunneling Microscopy 

1. Low Temperature Mode 

Nucleation of Ge clusters at low temperatures (< 600 °C) was recently described by us in 
a number of articles. ^'^'^^ We have demonstrated that huts form via parallel nucleation of 
two characteristic embryos different only in symmetry and composed by epitaxially oriented 
Ge dimer pairs and chains of four dimers on tops of the wetting layer M x N patches: an 
individual embryo for each species of huts — pyramids or wedges (Fig. [TK-c).'''^ These nuclei 
always arise on sufficiently large WL patches: there must be enough space for a nucleus on 
a single patch; a nucleus cannot form on more than one patch. This fact may be explained 
in assumption of presence of the Ehrlich-Schwoebel barriers^^ on sides of WL patches which 
prevent the spread of dimer chais composing nuclei from one patch to another; nuclei form 
on bottoms of potential wells of the spatial potential relief associated with the WL patches. 
As it follows from Ref.fiol. the total patch thickness (from Si/Ge interface to the patch top) 
rather than the mean thickness of WL controls the hut nucleation process. Further growth of 
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FIG. 1. (Color online) STM empty-state images of hut nuclei on Ge wetting layer formed at Tgj- — 
360°C: (a) pyramid (1) and wedge (2) nuclei on the adjacent M x N patches of wetting layer; hce = 
6 A; the structural models^*^ are superimposed on the corresponding image; (b) /iQe — 6A;^ (c) 
hce = 5.4 A. 



a hut likely decreases and finally eliminates the potential barrier; then hut occupies adjacent 
WL patches. 

The hut nuclei of both types arise at the same WL patch thickness of AML^-^^ at 
least at the initial phase of hut array formation, likely until growing huts start to contribute 
considerably to consumption of migrating adatoms.'' Therefore, they form at the same WL 
stress to relieve it. Consequently, they appear at the same strain energy (and with equal 
probabilities (P),-*^ hence, initially their free energy barriers (Vl^) are the same magnitude 
since these quantities are related through the Boltzmann factor:^^ P oc exp{—W/kTgj.}). 
This means that they are degenerate by the formation energy: if they had different formation 
energies they would first appear at different WL patch thicknesses.^^ Their nucleation rates 
are also equal at the initial phase of array formation {[dN^ / dt] / [dN^ / dt] = P^/P^^ w and p 
denote wedge and pyramid). 

Parallel nucleation of two competitive formations (hut nuclei) , different only in symmetry, 
on Ge WL is an amazing process. In Ref.liol. we proposed a speculation which somewhat 
explains this phenomenon on the basis of modeling of Ge cluster formation energy performed 



in Ref . 1221, according to which flat Ge islands likely occur on WL because of an energy benefit 
which arises in exposing the compressed {105} facets, rather than in relaxing the volumetric 
elastic energy. This speculation implies some difference in hut nucleation process from the 
usual the Stranski-Krastanow mechanism. However, this model certainly cannot explain 
the specific symmetries of nuclei and the fact that only two nuclei shapes appear on WL 



patches. 



10 the simultane- 



Another hypothesis might explain the latter facts. We reported in Ref. 
ous formation of the c(4 x 2) and p{2 x 2) reconstructions on tops of adjacent WL patches; 
a pyramid nucleus was observed on the surface of a c(4 x 2)-reconstructed patch. ^^^^^ We 
suppose that WL patch top reconstruction determines the type of a hat which can form 
on a patch. If this is the case, the c(4 x 2) reconstruction of the patch top enables the 
nucleation of a pyramid, whereas the p{2 x 2) reconstruction allows a wedge to arise on 
the patch. Domination of one of these reconstruction types may result in domination of a 
certain species of huts. 

Coexistence of the two reconstructions of WL patch tops seems to be more explicable. If 
they form simultaneously they have almost equal energies and correspond to two potential 
minima among possible variations of the (2x1) dimer-row structure (they have two simplest 
possible symmetries). Let us consider two adjacent dimer rows forming on top of a WL patch. 
A fluctuation may give rise to parallel orientation of buckling of two neighbouring dimers 
from different rows (when one of the closest atoms from different dimer rows is in upper 
position and another is in lower position) or result in opposite orientation of dimer buckling 
(when both closest atoms from different rows are in upper or in lower position). The next 
dimers in each row will stand crosswise to the flrst one to minimise the energy; they will form 
the inverse conflguration of pairs. Then the system should transit like a zipper to the state 
corresponding to one of two possible potential minima: p{2 x 2) for the case of the parallel 
buckling of initial dimers or c(4 x 2) for the opposite one. The resultant reconstruction will 
fill an entire top layer of the patch but it will not affect a state of the adjacent patch and any 
of the two possible reconstructions can arise on its top, that is observed experimentally.— 
These reconstructions may be stable at low temperatures of Ge deposition and each species 
of huts may nucleate on a patch with a specific structure of its top. 

Thus the, according to our assumption, fiuctuations of relative buckling of only two 
dimers in adjacent rows on top of each patch, which seem to be very likely, giving rise 
to simultaneous parallel formation of two modifications of the (2x1) structures with the 
simplest symmetries, may be the reason of parallel nucleation of the two known species of 
huts at low temperatures. 

Remark that this hypothesis certainly needs a theoretical and experimental verification. 
A usual model of the Stranski-Krastanow hut cluster formation will also require modification 
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if the hypothesis is true. 

In conclusion of this section, we would like to consider an alternative scenario of serial 
nucleation of huts which also may take place. Let us address again to Fig. [1^. One can see 
that the formation marked by the numeral 2 and considered as a wedge nucleus is composed 
by two couples of dimer rows, each row consists of four dimers (an arrow which starts at the 
numeral 2 indicates the lower right couple of rows). Such couples of 4-dimer sections of rows 
are often seen across the WL surface in the STM images .-^^i^^ They may be considered as 
cluster embryos in the process of serial nucleation. Such embryo may compose a nucleus of a 
wedge by pairing with another embryo and forming a structure marked by the numeral 2 in 
Fig. dK. Then, after reconstruction, the latter structure may transform in the blossom-like 
formation considered as a nucleus of a pyramid (marked by 1 in Fig. [TK). Perhaps, this 
scenario gives a direction of thought to explain the decrease in the fraction of pyramids in 
the arrays which is observed as the Ge coverage increases .-^^ 

We note also that this scenario does not contradict the above hypothesis about the affect 
of the patch top reconstruction on hut nucleation process; both reasonings can be easily 
combined to explain the whole process of hut nucleation. 

2. High Temperature Mode 

In contrast to the low-temperature mode, there are at least two ways of nucleation of 
{10 5} -faceted Ge clusters at high temperatures.— 

The first way is similar to the process of hut nucleation at low temperatures. Pyramids 
were observed to nucleate in such a way. Fig. [2K,b illustrates this statement: the pyramid 
nuclei, absolutely the same as those observed in the samples grown at low temperature, are 
seen on the WL patches in the images of the samples obtained at Tgr = 600 and 650°C, 
Hqq = 5 a. Their density was small, and they were mainly situated in the vicinity to large 
mature pyramids (Fig. [2]3,d), which arise at early stages of Ge deposition and have much 
greater sizes (up to 8 nm in height at /iQe = 6A)^^ than huts formed at low temperatures 
at the same values of hce-^^^ The WL surface mainly consisted of monoatomic steps and 
narrow terraces in these ares (Fig. [Jh). Neither wedges nor wedge nuclei were observed. 

The absence of wedges or even wedge nuclei at high temperatures, in contrast to low tem- 
peratures, waits for a theoretical explanation. It does not contradict the proposed qualitative 
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FIG. 2. (Color online) STM empty-state images of hut nuclei (a, b) and pyramids (c, d) on Ge 
wetting layer formed at high temperatures: (a, b) nuclei, hce = 5 A; (a) Tgr = 600°C, 43 x 37 nm, 
(1) indicates the pyramid nucleus; (b) Tgr = 650°C, 7.8 x 6nm; (c, d) pyramids, hce = 6A; (c) 
Tgr = 600°C, 110 X 110 nm; (d) Tgr = 650°C, 200 x 200 nm. 

speculations, however. 

The second way, somewhat resembling the process described by Goldfarb et al.^ for the 
case of the gas-source MBE (and thick hydrogenated WL), is illustrated by Fig. [3l At small 
values of /iQe, regions containing excess of Ge atoms are observed on the surface. Usually, 
they are not resolved as structured formations and resemble shapeless heaps of Ge (Fig.[3K). 
Terrace edges or pits^^ usually accompany them. Heap density is about lO^cm"^. Some 
of heaps start to form the {105} facets during Ge deposition (Fig. [3b). Annealing at Tgr 
results in formation of volumetric structures partially faceted by {105} planes, transforming 
heaps to some similarities of huts but with truncated very disordered and rough apexes 
(Fig. Eb-e). Formations of this type are also observed in Fig. [2l3,d: they are adjacent to 
mature pyramids. 

We have never observed such formations and such process at low temperatures of growth 
and suppose them to be inherent only to the high-temperature assembly of huts. We assume 

8 



FIG. 3. (Color online) Formation of {105} facets on shapeless heaps of Ge: STM images of different 
phases of faceting, Tgr = 650°C, hce = 5 A; (a) a shapeless Ge heap without faceting, 150 x 141 nm; 
(b) the outset of faceting, 64 x 64 nm; (c)-(e) effect of annealing at Tgr = 650°C: developed {105} 
facets are clearly seen; (c) 72 x 72 nm; (d) 46 x 46 nm; (e) 23 x 23 nm. 

that Ge heaps appearing at high temperatures are products of intense migration of Ge 
adatoms to surface irregularities, such as pits or terrace edges or sides of large pyramids, 
which are active attractors of highly mobile adatoms on Ge WL due to WL local stress 
relaxation. Much less mobility (the hop rate)^^ of adsorbed Ge atoms is likely a cause 
of absence of Ge heaps at low temperatures and formation of huts from the plane 16- 
dimer nuclei: huts nucleate before most of adatoms can reach sinks. Hut nucleation via 
formation of the 16-dimer embryos appears to be the most efficient sink of Ge adatoms at 
low temperatures, whereas at high temperatures, significant effect of competitive attractors 
redirect adatom fiuxes and decrease hut nucleation rate. 

An alternative explanation of the observed phenomenon may be as follows: Minor amount 
of carbon perhaps arriving to the surface and into the Ge film from the residual atmosphere of 
the MBE vessel or from chemicals during Si substrate preparation, or even from Czochralski- 
grown Si wafers might also give rise to formation of the Ge heaps via Ge adatoms migration 
to WL areas with local strain lowered due to presence of carbon. A possibility of presence of 
some difficult-to-detect amount of carbon in the Ge films grown at 600 or 650°C is hinted by 
formation of the c(4 x 4) reconstruction on the Si(OOl) surface during its dehydrogenation 
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FIG. 4. In situ RHEED patterns of Ge/Si(001) films, E = lOkeV, [110] azimuth: (a) Tg^ = 650°C, 
hce = 4A; (b) Tgr = 360°C, hce = 4A; (c) T^, = 650°C, hce = SA; (d) T^, = 650°C, hce = SA, 
annealing at the deposition temperature for 7 min; (e) Tgr = 650°C, /iCe = 6 A, the similar pattern 
is obtained for Tg^ = 600°C. 

at the temperatures higher than 600°C in the MBE chamber after wet deoxidization of 
the Si surface in dilute aqueous solutions of the hydrofluoric acid.^^ However, we have not 
got a veriflcation of this assumption at present and believe that it is very unlike. Further 
experiments on Ge fllm growth on buffer Si layers shall verify or refute this hypothesis. 

B. Reflected High-Energy Electron Diffraction 

Diffraction patterns of reflected high-energy electrons for samples of thin (/iQe = 4A) 
Ge/Si(001) fllms deposited at high (650 or 600°C) and low (360°C) temperatures with equal 
effective thicknesses are presented in Fig. [4K,b. The patterns are similar and represent 
a typical (2 x 1) structure of Ge WL; reflexes associated with appearance of huts (the 
3D-reflexes) are absent in both images, that agrees with the data of the STM analysis. 
Diffraction patterns presented in Fig. [4K,c,e correspond to the samples with hce increasing 
from 4 to 6 A. The 3D-reflexes are observed only in the pattern of the samples with hce = 
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6 A , that is also in good agreement with the STM data. 

Influence of the sample annealing at the deposition temperature is illustrated by a com- 
plimentary pair of the RHEED patterns given in Fig. [4l3,d. Annealing of specimens at the 
temperature of growth resulted in appearance of the 3D-reflexes (Fig. [4]1) that also corre- 
sponds to the results of our STM study. 

Difference in dynamics of diffraction patterns during the deposition of Ge atoms is a 
characteristic feature of the high-temperature mode of growth in comparison with the low- 
temperature regime. The original Si(OOl) surface is (2 x 1) reconstructed. At high temper- 
atures, as hce increases, diffraction patterns evolve as (2 x 1) ^ (1 x 1) ^ (2 x 1) with 
very weak ^-reflexes. Brightness of the ^-reflexes increases (the (2 x 1) structure becomes 
pronounced) and the 3D-reflexes arise only during sample cooling. At low temperatures, the 
structure changes as (2 x 1) ^ (1 x 1) ^ (2 x 1) ^ (2 x 1) + 3D-reflexes. This observation 
reflects the process of Ge cluster 'condensation' from the 2D gas of mobile Ge adatoms. 
High Ge mobility and low cluster nucleation rate in comparison with fluxes to competitive 
sinks of adatoms determines the observed difference in the surface structure formation at 
high temperatures as compared with that at low temperatures. 

IV. CONCLUSION 

Summarizing the above we can conclude that, as distinct from the case of low-temperature 
deposition of Ge fllms on Si(OOl) when the only process of Ge hut appearance exists, two 
process resulting in emergence of {105}-faceted clusters on Ge WL have been observed at 
high temperatures: Pyramids have been observed to nucleate via the previously described 
formation of strictly determined structures, resembling blossoms, composed by 16 dimers 
grouped in pairs and chains of 4 dimes on tops of the wetting layer M x N patches, each 
on top of a separate single patch, just like it goes on at low temperatures; an alternative 
process consists in faceting of shapeless heaps of excess Ge atoms which arise in the vicinity 
of strong sinks of adatoms, such as pits or steps. The latter process has never been observed 
at low temperatures; it is typical only for the high-temperature deposition mode. 

Different dynamics of RHEED patterns during the deposition of Ge atoms in different 
growth modes is observed, which reflects the difference in mobilities and adatom 'conden- 
sation' fluxes from Ge 2D gas on the surface for different modes, which in turn control the 
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nucleation rates and densities of Ge clusters. 
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